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Abstract 

This paper discusses the capabilities of the Interaction of 
Structures, Aerodynamics, and Controls (ISAC) system of 
program modules. The major modeling, analysis, and data 
management components of ISAC are identified. Equations of 
motion are displayed for a Laplace-domain representation of 
the unsteady aerodynamic forces. Options for approximating a 
frequency-domain (ico) representation of unsteady 
aerodynamic forces with rational functions of the Laplace 
variable (s) are shown. Linear time invariant state-space 
equations of motion that result are discussed. Model 
generation and analyses of stability and dynamic response 
characteristics are shown for an aeroelastic vehicle which 
illustrate some of the capabilities of ISAC as a modeling and 
analysis tool for aeroelastic applications. 

Introduction 

The goal of developing more profitable transport aircraft 
tends to result in more flexible and less inherently stable 
vehicles. The additional flexibility increases the likelihood of 
the need for active flutter suppression and gust and maneuver 
load alleviation. The lessened inherent stability results in the 
need for higher gain and, therefore, higher bandwidth ’’rigid 
mode" control. Both effects reduce the frequency separation 
between rigid and elastic modes. The desire for less 
observable, more agile fighters possessing good maneuvering 
capability over a large range of flight conditions leads to a 
need for similar analysis capabilities. Consequently, it is 
essential to have tools that consider the interactions that may 
arise between flexible structures, unsteady aerodynamics, and 
active controls. 

During the early 1970’s, the NASA Langley Research 
Center did not have access to an analysis tool to investigate 
such interactions. The Interaction of Structures, 
Aerodynamics, and Controls (ISAC) system of program 
modules was developed during the 1970's to provide the 
missing analysis capability 1 and continues to be an effective 
tool. Enhancements and refinements have been made during 
the 1980's. 2 * 3 * 4 * 5 Other aeroelastic packages containing 
similar analysis capabilities have been described in the 
literature. They include FLEXSTAB/* NASTRAN,? 
DYLOFLEX, 8 ADAM, 9 STARS, 10 and several industry- 
developed, proprietary codes. 

This paper describes and illustrates the application of the 
ISAC system of program modules. ISAC provides an 
integrated assembly of linear techniques for modeling, data 
management, and analysis of actively controlled aeroelastic 
vehicles. Typical applications of ISAC include: generation of 
sensor output equations and rational function approximations 
of unsteady aerodynamic data for use in support of simulations 
of aircraft; development of reduced order design plant models 
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for use in control law synthesis; evaluation of the performance 
of candidate control laws; as well as evaluation of the 
aeroelastic behavior of the open-loop vehicle. 

ISAC is operable in either a batch or an interactive mode 
and provides a variety of selectable graphical outputs with 
which to track model generation and assess dynamic behavior. 
Results are shown that were obtained by using ISAC to model 
and analyze a flexible, research vehicle. The results herein 
and those of the references document many of the modeling 
and analysis options provided by ISAC. 

Equations of Motion 

Figure 1 illustrates the fact that ISAC expresses aircraft 
motion as perturbations from steady -state conditions in terms 
of a superposition of characteristic mode shapes (including 
rigid body modes when appropriate). Bold symbols in the 
figures and text refer to matrix or vector quantities. Variables 
shown in figure 1 are the generalized coordinates (£), 
displacement at a point (x, y) of a modeshape (Zi), and 
displacement (z) at a point (x, y) due to a superposition of 
modeshapes. The modal characteristics, planform geometry, 
and flight conditions of interest are the basic inputs required 
by ISAC to initiate the development of aeroelastic models. 


z ( x, y, t ) = J) Zj ( x, y ) * ( t ) 

_ J l 

DISPLACEMENT GENERALIZED COORDINATES 

Fig. 1 Displacement as superposition of modeshapes. 

The primary dynamic blocks that are included in ISAC in 
the development of a closed-loop mathematical model are 
shown in figure 2. They are: flexible aircraft (an externally 
provided modal characterization of the flexible aircraft that is 
typically, but not necessarily, in terms of an orthogonal set of 
in vacuo modes); aerodynamics (generalized unsteady 
aerodynamic forces which are Mach number and frequency 
dependent); sensor, controller; actuator; and gust spectra 
dynamic characteristics. The variable, r\ , represents 
statistically independent white noise inputs which, when 
passed through the gust filter, produce gust incidence angles, 
<Xg, whose output power spectral densities correspond with 
user-defined gust spectra. Other, not previously defined, 
variables in figure 2 are the control surface deflections (5), the 
commands to the actuators (u), an external vector of 
commands (Ucom), and a vector of outputs (y). 






( c ) 8tj( s ) - Ssi ^7. the usual viscous damping form 
when the stiffness matrix is diagonal 
(m Wi®nj = k Wi md 8 ti (s) = 0 if ©„( = 0). 

The variable, g S j, is a structural damping coefficient. Only the 
third form is allowed in the state-space formulation below. 

Tg(s) is a diagonal matrix of actuator transfer function 
dynamics, and q is dynamic pressure. Ts(s), the sensor 

dynamics matrix, is also diagonal. The elements of C( s) in 
equation (2) are the contributions to the output y per unit 
generalized coordinate displacements, control deflections, and 
incremental incidence angles due to turbulence. Coordinate 
system options which admit analyses of perturbations from 
nonrectilinear flight conditions are also selectable in ISAC. 

Closed-loop analyses are accomplished by replacing u with 
the following equation 


Fig. 2 Dynamic elements included in equations of motion. 


U = -Tc(s) Pcyy + Ucom (3) 


Lap lace- Domain Form 

The Laplace-domain form of the plant equations of motion 
that are implemented in ISAC for an inertial, flutter-axis 
coordinate system and rectilinear motion are shown in 
equations (1) and (2). 


G^(s) G^s(s) 
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G^=M^s 2 -qQ^ , 


and 


y = Ts(s> C(s) V 

= T S (S) [C( (S) § + Cs (S) 5 + Cg(s) Og], (2) 

Here C(s) and are defined in terms of their components, s 
refers to the Laplace variable, and the subscript g refers to gust 
inputs. The variable Q = Q(Nma» sb/V) refers to generalized 
aerodynamic forces arising from motion of the generalized 
coordinates, control surface deflections, and gust velocity 
inputs. The variable, Nma* is Mach number, V is free stream 
airspeed, and b is a reference length, typically the mean 
aerodynamic semichord. The generalized mass (M^, M^&), 
generalized stiffness, K^, and modeshape information, 
together with the model for structural damping defined below, 
define the modal characterization of the structure. The symbol 
ktfr refers to the i* diagonal element of K». Up to three 
choices can be made by the user for the variable g^(s): 

(a) g^.(s) = V^Tgiit the usual structural damping 

(b) ^.(s) = jjjgji* reduces to (a) for s = ico, or 


where Pc y is a rectangular matrix of ones and zeros that 
selects a subset of the available outputs for feedback, Tc(s) is 
a transfer matrix representation of a multi- input/multi-output 
controller, and u CO m is a vector of external commands to the 
actuators. 

The discussion thus far has assumed the availability of 
unsteady aerodynamic data in the complex s-plane. The 
aerodynamic data is typically only available on the imaginary 
axis at a finite number of frequencies. Nevertheless, with the 
Laplace-domain form of the equations and the approximation 
that Q(s) * Q(0 + itob/V) for lightly damped systems, 
commonly referred to as the p-k approximation, 11 stability 
analyses for open- and closed-loop systems can be performed 
using determinant or matrix iteration.^* 2 * 1 1 Frequency -domain 
analyses can also be performed to exhibit frequency responses 
and (for stable systems) power spectral densities and rms 
results for outputs of interest. Time-domain analyses can be 
performed using fast Fourier transform (FFT) techniques. 12 
Results illustrating these capabilities will be shown below. 

Linear Time Invariant State-Space Form 

ISAC also provides the capability of converting the 
equations to a first order linear time invariant (LTI) state-space 
form with subsequent stability analyses performed using 
efficient, linear algebra techniques. This conversion is 
accomplished by approximating the frequency-dependent 
unsteady aerodynamic forces that are defined for a set of 
points on the s-plane imaginary axis with rational functions in 
s, the Laplace variable.^ 13,14,15,16, 17 Two rational function 
approximation options available within ISAC to use in making 
the conversion to LTI state-space equations are shown in 
equations (4) and (5): 


6-Ao+ A,^ +A 2 (y) + X A 2+j — ^vai W 

M s+ b 

or, for R diagonal, positive definite, and x ^2. 

Q=Ao+ A,* +A 2 (y) 2 + d[yI - rJ * E. (5) 
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In these equations, the matrix approximate, Q, is 
n£X(n^ + n$ + ng); and l\ and I 2 are the number of unique 
denominator poles in each type of approximation, respectively. 

The Ai, 6, R and € are real, constant matrices, and each 
(-Vpj fb) is a stable, real denominator pole. The matrices D 

and E are of rank / 2 for the usual case where n^ > £ 2 - 

The Roger’s 13 form, equation (4), allows rapid solution for 
the coefficient matrices that best fit the tabular data in a least 
squares sense, particularly if the denominator coefficients are 
chosen a priori rather than by optimization of a performance 
measure. However, the dimension of the added states 
corresponding to the rational approximation can be large, 
being equal to the product (n^ x / 1 ). 

In the minimum-state approach 16 (equation (5)), because 

the matrices D (n^ x l 2 ) and E {l 2 x (n^ + n 6 + ng» are of rank 
I 2 , only 1 2 added states are required, independent of the 
number of generalized coordinates. Therefore, the minimum - 
state approach can offer a significant savings in number of 
added states with little or no penalty in accuracy of modeling 
the aerodynamic forces. However, solution for the coefficient 
matrices is time intensive requiring an iterative, two-step 
process. 

A third option, not explicitly presented here but available in 
IS AC, allows a columnwise fit 5 * 17 of the form of equation (4). 
Each column can then have a different number of lag terms 
with different values for the denominator poles. A discussion 
of the advantages and disadvantages of these options is given 
in references. Reference 15 discusses an approach, not 
currently implemented in ISAC, which, given the 
characteristic poles of the system for a specific flight 
condition, allows the rational portion of the unsteady 
aerodynamic forces to be modeled with no added state penalty. 

These approximations, together with rational 
representations of sensor, actuator and gust spectra dynamics, 
allow a state-space model of the aeroelastic vehicle to be 
generated of the form shown in equations (6) and (7). 


and 


x = Ax + B u u + Bgt], (6) 

y =Cx + D u u+ DgT|. (7) 

See the appendix for more detail. The state vector is 


* T =(f.4 T .«*A 


X>- 


The subscripts a, S, A, and g refer to states due the rational 
function approximations, sensor dynamics, actuator dynamics 
and the gust filter dynamics, respectively. 

The D u and Dg can arise because of the rational function 
approximations. Consider the case of an accelerometer with 
negligible dynamics in the frequency range of interest. If 
s 2 Tsj is strictly proper (that is, the numerator is of lessor 
degree than the denominator) for each of the actuators, then D u 
will be zero. If this condition is violated for the jth actuator, a 
nonzero D u can be avoided by constraining to zero that 
column of A 2 (see eqs. (4) and (5)) which would make s^Tgj 
not strictly proper. Once a choice is made for gust spectra 
dynamics, similar considerations would be necessary to avoid 
the presence of Dg in an accelerometer output. 

Once the state- space form is generated, it, or a reduced 
order approximation thereof, can be used as a design model of 
the plant to perform controller synthesis studies using state- 
space design techniques. Alternatively, tabular frequency - 
response representations developed without the need for the 


s-plane approximations can be used as design models for 
frequency-domain-based design techniques. 18 

Closed-loop analyses with the state-space representation 
can be accomplished by direct state- space definition of the 
controller or by converting a transfer-matrix controller 
representation to a state-space equivalent. The direct state- 
space controller is assumed to be of the form 


XC = AcXc + BcPcyy. 

(8) 

U = (CcXc + D c Pc y y)+ Ucom. 

(9) 


where the subscript C refers to the controller. The matrix Pcy 
selects a subset of the available outputs to use as feedback. 

Mqior ISAC Components 

The major ISAC components and their data communication 
links via the data complex are identified in figure 3; the figure 
also shows how ISAC outputs can provide aeroelastic design 
models for control law synthesis. 



Fig. 3 Major modules contained in ISAC. 

Data Storage and Mana£ement 

The model development process is initiated by receipt and 
storage on the data complex of modal data from an external 
source. The data complex manager (DCM) facilitates this task 
and queries the user for dimensional and format characteristics 
of the data. The DCM catalogs the data with a user-provided 
descriptor and adds dimensional descriptors together with the 
time of storage. The DCM and the data complex are vital 
parts of the ISAC system of modules. The data complex is a 
random access storage facility for binary data. The data stored 
can include thirteen different types of data that represent not 
only information required to perform the aeroelastic modeling 
but also the state-space quadruples for each model that results. 
A single data complex may contain up to 100 separate data 
sets, each of which may contain all thirteen data types. Inputs 
made by the user direct the various ISAC modules as to where 
to store data received/created and how to uniquely identify that 
data for subsequent use. In addition, the DCM provides a user 
friendly interface allowing direct interactive access to the data 
complex contents. 

Unsteady Aerodynamic Force Computations 

The two modules DLIN and DLAT generate the basic 
information required to include unsteady aerodynamic forces 
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in the mathematical models. DLIN is a modal interpolation 
preprocessor for the DLAT module; and DLAT calculates 
unsteady aerodynamic forces using a doublet lattice 19 method. 
Alternatively, one or both of these modules could be bypassed 
if the data they generate were available from an external 
source. Once the aerodynamic description is complete, model 
generation at flight conditions of interest can be accomplished 
to provide design models or to perform analyses of stability 
and dynamic response characteristics. 

The DLIN module accepts modeshape data from the data 
complex or from an external source together with information 
related to that data such as the number of structural sections 
and whether each structural section corresponds to a plate or a 
beam representation. A spline fit is made of the modeshape 
data. The next step is to utilize aerodynamic paneling inputs 
to determine the points at which interpolation is needed to 
obtain modal deflections and/or slopes. The interpolations are 
performed and the resulting data are stored on the data 
complex for subsequent use by the DLAT module in 
computing generalized aerodynamic forces. 

The DLAT module is employed to compute unsteady 
aerodynamic force data for a specified Mach number and for a 
specified set of reduced (nondimensional) frequencies. The 
aerodynamic force data can be computed and saved on the data 
complex in two forms, Q aic(Nma, icob/V) and 
Q(Nma» itob/V). The matrix of modeshape independent 
aerodynamic influence coefficients, Qaic (Nma. itob/V), is of 
dimension nb x nb, where nb is the number of boxes in the 
aerodynamic paneling. The ijth element of the aerodynamic 
influence matrix contains the pressure at box i due to a normal 
component of motion of box j. The matrix of generalized 
aerodynamic forces, Q(N M a. icob/V), is of dimension 
n$ x (n$ + n$ + ng). The aerodynamic influence coefficient 
matrix, together with modal data (e. g., that provided by the 
DLIN module), can be used to compute generalized forces as 
shown in equation (10). 

Qij(NMA.^) = ZTQa3c(Nma, ^) wj, (10) 

where Qh is the contribution to the generalized aerodynamic 
force in the ith degree of freedom due to harmonic motion in 
the jth degree of freedom, Zi is displacement in the ith mode at 
each box pressure point, and wj is down wash at each box 
downwash point due to harmonic motion of the jth mode. 

A selectable set of rigid body, gust, and control rotation 
modes are also available within the DLAT module and can be 
employed to obtain stability and control derivative information 
and/or to augment elastic mode information from other 
sources. The doublet lattice 19 procedure is an integral part of 
the DLAT module (Augmentation is planned to enable 
generation of piston theory aerodynamic forces within ISAC 
for high Mach number unsteady conditions). 

Model Formation and Analysis 

Components used to determine stability and response 
characteristics as well as to generate s-plane approximations of 
the unsteady aerodynamic forces are contained within a 
module known as DYNARES. The primary functions of 
DYNARES are to assemble a model using the data previously 
described, to perform open- or closed-loop analyses of 
dynamic characteristics, and to provide reduced order plant 
design models for use in control law synthesis. Analyses 
performed with the model in the form of equations (1) and (2) 
can use the tabular aerodynamic data from DLAT or, 
alternatively, the rational function approximations thereof. 
Stability analyses using equations (1) and (2) are performed by 
the module STABCAR,^ which is not directly connected to the 


data complex; but the DCM module can be used to provide an 
input data file for STABCAR. 

The s-plane capabilities within DYNARES can be used to 
determine the rational function approximations to the 
frequency -dependent aerodynamic data using previously 
discussed Roger's or minimum -state approaches provided the 
denominator poles are chosen a priori. The data representing 
the resulting rational function approximations are stored on the 
data complex for use in forming a LTI state-space 
representation. Alternatively, a separate module, SPLFIT, can 
be employed to determine the column- independent, column- 
dependent, or minimum-state rational function 
approximations. This module also provides for optimization 
of the denominator pole locations. ^ The resulting data can be 
stored onto the data complex for subsequent use. 

Sample Application Results 

Results obtained with ISAC will now be presented that 
illustrate its utility in mathematical model development and in 
aeroelastic analyses. 

Example Vehicle Description 

The vehicle modeled and analyzed is an unpiloted drone, 
the DAST ARW-2 (drones for aerodynamic and structural 
testing, aeroelastic research wing number two) 20 * 21 shown in 
figure 4. All-moving stabilizers are identified and the 
locations of inboard and outboard ailerons are also indicated. 
A subset of the vehicle sensors and their locations are noted on 
the figure. The wing design and placement on the fuselage 
was such that flutter suppression, load alleviation, and rigid 
mode stability augmentation were required in regions of the 
vehicle flight envelope. 
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Fig. 5 Paneling and structural node information. 



Inputs for Ae rodynamic Force Computations 

Figure 5 shows the planform views of the wing, horizontal 
tail, vertical tail, and a portion of the fuselage modeled as a 
lifting body. Alternatively, the fuselage could have been 
modeled as a slender body. For convenience, the vertical tail 
is shown in this view by translating it outboard and displaying 
its side view. 

The boxes depict the elemental areas required for a doublet 
lattice representation of the lifting surfaces in terms of 
doublets on the one-quarter chord lines of each box. The 
boundaries of the control surfaces coincide with elemental box 
edges and are shown in the figure as darkened lines. There are 
no darkened lines on the vertical tail because it has no rudder. 

The circles in figure 5 denote points at which modal data 
were available and indicate that the wing was modeled as a 
plate; whereas the fuselage and tails were modeled as beams. 
Twelve symmetric degrees of freedom were considered in the 
results to follow (plunge, pitch, and ten elastic modes) together 
with three control modes and one vertical gust input. At the 
mid span of each elemental box, displacement (at the one- 
quarter chord point) and displacement and slope (at the three- 
quarter chord point) are required for each generalized 
coordinate, control deflection and gust incidence angle in 
order to generate generalized aerodynamic forces 
corresponding to a particular Mach number and 
nondimensional frequency. Cubic splines in one (for beams) 
and two (for plates) dimensions 22 are used to perform the 
interpolations for modal deflections and slopes at the points 
required by the doublet lattice code. Examination of figure 5 
and graphical display of modal deflection data, such as shown 
in figure 1, provide a valuable check for the correctness of the 
inputs into the DLAT module that are required for 
aerodynamic force generation prior to initiation of the 
computations. 

Rational Fun ction Approximations 

Approximations of two elements of the unsteady 
aerodynamic force matrices are shown in figure 6 that were 
obtained using the Roger’s option with l\ - 2. The elements 
chosen are from the gust column, typically the column most 
difficult to fit. The figure provides a polar depiction of how 
the elements vary as a function of the (hidden) frequency 
variation; it compares the tabular data and the corresponding 
rational function approximations. The direction of increasing 
frequency is indicated on the figure. The approximation was 
constrained to fit the data precisely at zero frequency and to fit 



Fig. 6 Polar plots of generalized aero forces due to gust. 



b) Qgag 


Fig. 6 Concluded. 

the rate of change with respect to frequency at zero frequency. 
References 3 and 5 discuss the types of constraints that may be 
imposed in ISAC to ensure a good fit of the data in critical 
frequency regions. 

Stability and Response Characteristics 

Figure 7 shows the loci of the critical subset of elastic roots 
as a function of normalized controller feedback gain (the 
nominal gain has a magnitude of one) for a candidate single- 
input/single-output flutter suppression control law. 
Examination of the locus for the 4th mode reveals that the law 
does not achieve -6dB of gain margin. There are two curves 
for each root locus. These curves contrast the root loci 
obtained using the p-k approximation with root loci obtained 
using the rational function approximations for the unsteady 
aerodynamic forces. When these curves are generated using 
the iterative, Laplace-domain approach, 3 one can easily add a 
term in the controller transfer function corresponding to phase 
errors thereby exhibiting whether phase margin constraints are 
satisfied. For example, replace Tc(s) with exp(i7i/4) Tc(s) and 
repeat the analysis; if the system is stable at nominal gain, one 
has at least 45 degrees of margin to uniform phase errors in all 
inputs or in all outputs. Automated generation and graphical 
display of root loci are selectable within ISAC for feedback 
gain, density, altitude, and velocity variations under the 
assumption that Mach number is fixed. 



Fig. 7 Feedback gain loci of critical elastic roots. 
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(13) 



An illustration of gust response analyses performed using 
the frequency-domain form of the equations is shown in 
figure 8. This analysis, valid only when the system (open- or 
closed- loop) is stable, shows the power spectral density (psd) 
of the incremental wing root torsion moment due to turbulence 
having the Von K arm an spectrum. The psd shown has been 
normalized by lower dynamic pressure results. Also shown on 
the same figure is the cumulative rms of this load with 
frequency. These types of information can also be generated 
for other selectable outputs including control deflections and 
rates for Von Karman, Dry den, or user-defined rational 
polynomial definitions of turbulence spectra. 

Reduced order approximations of mathematical models are 
often desirable in order to isolate the critical interactions, 
provide smaller, more manageable state-space models for 
control law design, or reduce the time required for analysis. 
Within IS AC, truncation and residualization options are 
available and are applied in the Laplace-domain representation 
(see equations (1) and (2)) as follows: Use the second block of 
rows of equation (1) to eliminate 8 in the first block of rows; 
then partition § into low frequency, $ 1 . and stable, high 
frequency, £ 2 , modes yielding 


Gn(s) 
- G2i(s) 


G«(S) 1 


G22(S) . 1 

M) 

/ ^ \ 

[ qQ?g(s) (qQs8(s)-lVV2)T 8 (s) ] 8 

' F,(s) ‘ 

(+) 

- F 2 (s) . 

• r 


(ID 


Now retain only the static effects of £2 by assuming £2 reaches 
steady state instantaneously relative to £ 1 . Then 


I" G„(s) Gi 2 (0) ‘ 

5i 


' Fi(s) ‘ 


L G2!(S) 022(0) . 

^2, 


- F2(S) . 

UJ 


( 12 ) 


can be used to eliminate £ 2 - The elimination is accomplished 

by solving for £2 in terms of the other variables in the second 
block of rows of equation (12), substituting the result in the 
first block of rows, and collecting like terms. The result has 
the form 


G,i(s)i, =Fi(“ 8 ) 


and, by making the substitution for £2 in equation (2), one can 
obtain 

y = Ts(s) [CiGOl, + Cs(s)T6(s)u + Cg(s)Og]. (14) 

If LTT state-space models (eqs. (6) and (7)) are created after 
residualization, model outputs are generated using eq. (2) with 

Cl truncated, rather than eq. (14), which is an approximation 

if Ci(s = 0) * 0. Approximation can be avoided by 
performing residualization outside ISAC after generating a 
full-order state-space model. 

Results are shown in figure 9 from a combined truncation 
and residualization model reduction at a flight condition well 
beyond flutter that was used as the primary active flutter 
suppression design point. Rigid body and essentially 
noninteracting elastic modes below the frequency of 
flutter (120 rad/sec) were removed by truncation. Modes of 
higher frequency that only had small impact upon the flutter 
mechanism were residualized. This figure exhibits magnitude 
of the frequency response of the difference between two 
accelerometer outputs due to a commanded trailing edge 
outboard control surface input (Ucomieo) for full twelve 
mode model and for the five mode approximation. Good 
agreement is evident over a wide frequency range including 
the flutter frequency. 

Figure 9 also serves to illustrate the capability for 
generation and display of frequency responses of interest. 
These frequency responses can then be used in design of 
control laws. When experimentally derived frequency 
responses are available over a portion of the desired frequency 
range, one can combine the experimental/analytical data to 
improve the representation of the plant. 18 

Time responses are computed within ISAC using the 
frequency-domain form of the equations coupled with Fourier 
transform techniques. An outline of the approach follows for a 
commanded input: Compute the fast Fourier transform (FFT), 
denoted by (F in equation (15), of the input for a sample rate 
high enough to keep digital effects small in the frequency 
range of interest and to preclude aliasing within the user- 
controlled input frequency content. 

Ucomj(<0)= (Ucomj(t)). (15) 
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The frequency and time variables are shown as bold to 
emphasize that there are a set of values of Ucom^co) and 
Ucomjtf). Using equations (1), (2), and (3), compute the 
closed-loop frequency response, Hj(io>), of the outputs of 
interest. Hj(ico) need only be computed at a sufficient number 
of frequencies to allow accurate interpolation at frequencies 
corresponding to the FFT evaluation points where Uco^cd) is 
nonzero. Form the product 

y(a>) = Hj(ia>) Ucom/©))- (16) 

Compute the inverse FFT of the outputs. 

y(t)= <F- l (y(©)). (17) 

The time history generation process is illustrated in 
figures 10, 11, and 12 for symmetric input of a one-degree 
magnitude swept sine (sine chirp) command into each of the 
pair of trailing edge outboard control surface actuators 
(Ucomieo)' Figure 10 shows the computed frequency response 
over the frequency range of the input. The commanded input 
is shown in figure 1 1 where the input frequency varies 
logarithmically with time from 5 Hz to 20 Hz over a 
14.8 second time period. The amplitude of the input is ramped 
from zero to one degree maximum amplitude over a 
0.05 second interval, ramped out near the end over the same 
interval of time, and followed with a short period of zero 
input. The output time history that results from the control 
input is shown in figure 12. The results are for an open-loop 
condition well below the flutter boundary. The composite 
sensor output is the difference of two outboard accelerometers 
located at approximately the same spanwise location (see 
figure 4). Consequently, the composite sensor output 
accentuates the observability of torsional motion. 
"Blossoming” of the accelerometer time response is evident as 
the input frequency passes through the frequency of one of the 
lightly damped poles that form the flutter mechanism at a 
higher frequency and at a higher dynamic pressure. Response 
is also seen at higher frequencies corresponding to other lowly 
damped elastic modes. 




i 1 1 — * — * r 



- -02 1 i i i 1 I i i I I I 

0 5 in 15 

time, sec 

Fig. 1 1 Sine chirp command into outboard aileron actuator. 



The majority of the graphical outputs that have been 
presented pertain to the frequency -domain representation of 
the model. The comparison of the s-plane approximation and 
the underlying numerically generated frequency -domain 
aerodynamic forces (figure 6) and the automated generation of 
the loci of roots with gain (figure 7), density, altitude, or 
velocity variation were examples of state-space related 
analysis capabilities within ISAC. The state-space option is 
primarily used for stability analyses and in the generation of 
plant models for use in control law design by other tools such 
as MATLAB 23 and MATRIX*. 24 These external tools have 
capabilities that allow the generation of frequency and time 
responses of interest using the plant models. 

ISAC Portability. Availability, and Documentation 

ISAC is written in FORTRAN and will execute, with 
minimal changes relative to memory management, on many 
computer systems with a FORTRAN 77 compiler that supports 
namelist input. It is currently running on VAX micro II 






computers with VMS, SUN Sparc workstations with a UNIX 
operating system (4.1 or greater), and IBM minicomputers. It 
is available through NASA with its own plotting package and 
library routines, and with preliminary versions of a user’s 
manual and an extensive set of sample cases with 
documentation. 

The ISAC system of programs has been or is currently 
being used in support of numerous projects such as: 

- DAST ARW-125 and ARW-2^ 26 

- DC-10 wind tunnel flutter model 27 

- Generic X-wing feasibility studies 28 

- Analyses of elastic, oblique- wing aircraft 29 

- Active Flexible Wing (AFW) wind tunnel test 
program 1 8 *30.3 1 

- Generic hypersonic vehicles 32 * 33 

- Benchmark active controls testing project 

- High speed civil transport 

There is a growing user community. The ISAC system has 
been transmitted, upon request, to a number of research 
facilities which includes: 

- NASA Centers 

Langley Research Center 
Right Systems Directorate 

Guidance and Control Division 
Structures Directorate 

Structural Dynamics Division 
Dryden Flight Research Facility 

- United States Militaiy 

Department of the Air Force 

AFWAL at Wright Patterson AFB 
Department of the Navy 
NADC at Warminster 

- Companies 

Lockheed Engineering and Sciences Company 
Boeing Aircraft Company 
McDonnell Douglas Aircraft Company 
Learjet, Inc. 

Leth and Associates 

- Foreign Countries 

Australia 

Aeronautical Research Laboratories, Dept, 
of Defence 
Canada 

Aerodynamics/Dynamics of Right Directorate, 
Canadian National Defence Headquarters, 

Institute for Aerospace Research 

- Universities 

Technion University (Israel) 

Massachusetts Institute of Technology 
University of Washington, Seattle 
University of Rorida 

Concluding Remarks 

An overview has been presented of the theoretical bases, 
aeroservoelastic modeling and analysis capabilities, and 
facilities for data storage and management for the ISAC 
system of program modules. Results are presented for a 
flexible aircraft that illustrate the major ISAC capabilities. 
References are cited which present research performed using 
ISAC for a variety of configurations. The results presented 
here and in the references demonstrate that ISAC is a useful 
tool for aeroservoelastic modeling and analysis. 

ISAC is written in FORTRAN 77 and will execute, with 
minimal changes related to memory management, on many 
systems with a FORTRAN 77 compiler that supports namelist 
input. Preliminary versions of a user’s manual and 


documentation containing an extensive set of sample cases are 
available. 

Appendix; State-Space Equations in Flutter Axes 

The plant state space matrices for the flutter axes choice of 
coordinate system and Roger’s form (eq. (4)) of rational 
function approximations to the unsteady aerodynamic forces 
are more explicitly defined in this appendix. 
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= -5M- 1 * = -‘M- 1 (k,£ -qA§), 

A «~ 

See equation (4) for the Ai. The matrix is diagonal and 
viscous in form; the ith diagonal term is 

D 5i5i = S«i m 5&®n i - 

Refer to the expansion of terms below equation (1) for 
definition of the variables. 

1 1-1 0 0 “ 

0 11—1 0 — 


where each of the m rows contains i\ ones in the indicated 
columns. The aerodynamic states are ordered such that the 
effects of all lags upon each particular generalized coordinate 
are grouped together. 

A*—*" 1 !** °S *«] 
where 

^6 = -< lAo. Dj = -q-^Af , 0A.i = Mjij - A®. 

See the discussion of the actuator state rows below for a 
definition of and 



where the dynamic element whose output has the desired gust 
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velocity spectrum is required to have a denominator 
polynomial of order at least one greater than its numerator 

polynomial. The constraint A£ = 0 is also imposed. For the 
Dry den vortical gust spectrum 



where g L is a gust reference length and <x w is rms gust 
velocity. 






■j,— 

[ a «£ i a<l4 i i B °*] = 

'*« i (>A% | i('A«)C B | i(‘A*)D B 


; { * f 


where 


*A = 


ith row of A 3 
ith row of A 4 


ith row of A , 1+2 


are zero unless the sensors are on a control surface. The 
matrix with the g subscript allows inclusion of a gust vane 
sensor. 

The vector, y p , provides the inputs to the following sensor 
dynamics equations: 

*s = ^s x s + ®syp* ys = c s*s+ D sy P * 

The sensor system matrix is block diagonal with one block for 
each sensor. Currently, the system matrix for each sensor is 
expressed in numerically nonoptimum controller canonical 
form. The matrices in the sensor dynamics states rows are 

= +C/Sp)» A s £ = B 5 (CppA^ +C pp ), 

As* = BsC^A^, A ss = As , 

Asa =Bs(C>^ +C^C„ 2 +C* s C A[ +C D pS C A ), 

A*=B S (C>^ + C^). 

Bs« = Bs<C> 4m + CpgDfy ), E Sg = BsC> 6 - 

See the next paragraph for definition of D ^ 2 ■ 

The actuator state equations have block diagonal system 
matrices with each block representing one actuator in 
controller canonical form. It is required that each actuator 
transfer function have denominator polynomial at least two 
orders higher than its numerator. The following equations 
result: 


x a=AaXa + Ba“. » = CaXa, 

8 = Ca*a = C a A a x a = C^x^, 

8 =^!*^ =C A 1 ( A A x A + B A u ) = C A 2 x A +D A 2 u ' 

a aA = A a * B Au =B a . 


The gust state-space equations, defined previously, 
introduce the following matrices in the plant state-space 
equations: 



where there are n^ repetitions of the set of lag coefficients. 

The input into the sensors (output from the plant) has the 
form 


The following vector of output equations is generated in the 
state-space representation: 


y». “Ln>y£.( u - “com ) r . 


W&'iT, 



The vector L is a modal displacement characterization of 
incremental loads wherein 

L = c,£ 

The equations presented in the body of the text and in this 
appendix are sufficient to define the system matrices C, D u 
and Dg of equation (7). 


>P = C* 5 * C^S t C* $ t c',8 + c » + c “ 5 + 


The A refers to acceleration, the R to rate, and the D to 
displacement. The matrix coefficients with 5 in the subscript 
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